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ABSTRACT 


A computer code, FEASBL, is developed to maximize a 
non-linear objective function over a convex feasible region. 
The feasible region is defined by a set of non-linear and 
linear constraints on the variables.  FEASBL can solve prob- 
lems involving up to fifty variables with a feasible region 
formed by up to fifty non-linear constraints, and fifty 
linear constraints.  FEASBL uses a feasible direction method 


as its solution algorithm. 
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I. INTRODUCTION 


This thesis presents a computer code for solving medium 
sized non-linear programming problems, and describes its use 
at the Naval Postgraduate School computer center. 

The general non-linear programming problem is the optimi- 
zation of an objective function over a constrained feasible 
region. By the technique of separable programming, many non- 
linear problems can be converted to a linear problem and 
solved using established linear programming methods. Many 
problems, however, cannot be solved by this method. Several 
techniques have been developed for solving non-linear prob- 
lems. Most methods involve following the gradient of the 
objective function to generate a sequence of points that 
ΙΙ. το the Optimal point. The method used in this 
code is Zoutendijk's method of feasible directions [1]. The 
method of feasible directions can be considered a large-step 
gradient method as opposed to the small-step gradient and 
search methods.  Feasible direction methods can be charac- 
terized by: 

ШЕ сватешир ас ап initial feasible point. 

2. Establishing a feasible direction that will improve 
mm objective function. 

3. Proceeding to a new feasible point that has a better 
value of the objective function than the previous point. 

4. Converging to the optimal point, or determining that 


the problem is unbounded. 


Many procedures can be considered feasible directions methods, 
for example, the simplex method for solving linear program- 
ming problems. 

The method used in this code requires that the user know 
musobjective funebrono.cconstraint functions, all their first 
partial derivatives, and an initial feasible point. If the user 
does not know an initial feasible point, this code can be 
used to solve a "EE | programming problem which will 
provide a feasible point. This problem will be described 
later. This code will solve maximization problems of the 


form: 


Max: z = F(x) 


Subj to: G; (x) 


A 
o” 
Ro 
(n 
к= 

Pl 


aM οτε 
1 -- 


Enecrunctrions F(Xx);, and G. (x) аге поп-11пеат functions which 
are supplied by the user as Fortran FUNCTION routines. The 
sets Іс апа I, are the sets of subscripts for the non-linear, 
and linear constraints respectively. If the objective func- 
tion is concave, and the non-linear constraints are convex 
functions, then the problem will be a convex programming 
problem, and the method will converge to the global maximum. 
If the problem is not convex, then the method may converge 

to a local maximum. The code is designed to accommodate 


ο Шр to fifty variables. The constraint set 





can include up to fifty non-linear constraints and fifty 
linear constraints. A complete description of how to use 


this code is provided in Section IV. 





II. MATHEMATICAL BACKGROUND 


The theoretical foundation for the method of feasible 
directions was developed by Zoutendijk in 1960 [1]. A brief 
discussion of the mathematical background as it applies to 


the method used in this computer code is presented. 


A. DEFINITIONS 

υπ R iS a convex region if for Xi and Хо» two points 
in R, then Ehe point X4, = 8Xi > | T= а)х, is in R for all a 
such that 0 « a « 1. 

Let R be the feasible region, and if x is a member of R 
ШЕП а direction 5 is а feasible direction if and only if 
Ве хе а ^' » 0 such that x * As is in R for all 
ше XA'), 

Duc reassbiegdwerection is a feasible direction s 
such that v F(x)! s > 0. In other words, a usable feasible 
direction at x is a direction in which the value of the ob- 
jective function increases and along which.it is possible to 
move while remaining in R. 

ВЕ а concave function, and let the feasible 
region R be convex. Then if there exists no usable feasible 
ЗЕ СЕО ас х іп R; then x must be a maximum of F(x) over 
the region R. To show this assume that x is not a maximum 
um no usable feasible direction at x. There exists 
α роет зи В such that F(y) > F(x) and since R is convex 


NEMPE long the line between y and x must be in R. 
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Further, since F( ) is concave, v F(x) (y ΧΕ; 
emedi rection y - X is a usable feasible direction. This 


contradiction implies that x must be a maximum of F(x). 


Б ІНЕ METHOD 

The method of feasible directions begins at an initial 
feasible point хо, determines the 'best" usable feasible 
direction, moves along that direction to the point xÉ that 
maximizes F(x) along the direction without leaving the fea- 
Ге терісп, repeats this procedure until a point x* is 
found which has no usable direction, or determines that the 
problem is unbounded. The method can thus be viewed as a 


four step procedure. 


1. Determine a usable feasible direction SK at thc 
point хк. 
2. If no usable feasible direction exists stop, x“ 1$ 


the maximum of F(x) over R. 

SO etermine a step length Ак such that еті = хк + το; 
1s іп R and К maximizes F(x) along S 

фас ШЕ A. goes to infinity, then the problem is unbounded, 
Otherwise go back to step 1l and repeat. 

The direction finding problem will be discussed first, 
and then the step length problem. 

1. The Direction Finding Problem 

The direction finding problem is to find the usable 

feasible direction that will result in the greatest apparent 


increase in the objective function. This can be formulated 


as the following linear programming .problen. 
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Мах: σ 


μη το: o — v F(x) s Em 


Т 
g + V6, (X) 5 = 


A 
O 
“ 


iel (x) 


T 
а; $ < 0, lel, (x) 


5. > 0, jeJ (x) 


Бит, jeJ" (x) 


л 
^ 


Ех», 


мһете 1. (x) and I. (x) РЕВ Of binding non-linear con- 
straints and binding linear constraints, and J (x) and J” (x) 
Bre the sets of j such that x is at its lower and upper 
Pound respectively. 

Mene xs estat least one point that is strictly 
ui ror to R with respect to the non-linear constraints and 
1f the optimal value of c is O0 in the direction finding prob- 
lem, then there is no usable feasible direction and we are 
шие шенор ваши" розпе for the non-linear problem. The last 
constraint in the direction finding problem is a normalization 
Eustramt. In order to insure that the procedure will con- 
Bere the components of SK must be bounded, and several nor- 
malizations have been studied by Zoutendijk. These normalizations 


ате: 
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№3: ПТУ Ех). >0 
Er x 09, 


. 2-1 if У. Е(х). < 0 
i E], 


Т 
№: VG, (X) 5 % с < b; - G. (x); iel (x) 


Iw 


T | 
< b: га. х; iel, (x) 


The normalization Nl will generally lead to fewer steps then 
№, №, ог №, but the amount of computation for each step 
is considerably greater. Also the simple constraints в. < 0, 
and ze ο ο τν ο szıhessjzezo: the direction finding prob- 
lem if №1 is used, but not for N2, N3, and N4. As one goes 
from N2, to N4 the number of steps needed increases, but the 
amount of computation for each step decreases. N5 leads to 
a larger problem, and also does not converge as fast as the 
other normalizations. Normalization N2 was chosen for this 
code as the best compromise between the number of steps and 
the amount of work for each step. 
ο E ore calme Finding Problem 

To insure that the procedure actually converges to 
the maximum, some precautions must be taken to prevent the 
phenomenon of "jamming" in which the procedure in effect 
pets "stuck in a corner" and never finds the real maximum. 
This is accomplished by considering nearby constraints in 
addition to the binding constraints when solving the direc- 


tion finding problem. 
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Define the sets: 


Ic(x,t) = ЕТ С : b: - Е < G; (x) ep 
= (3 - T 
I; (x,€) - liel, : b: e<a,x < b.) 
Е Aye) t d. « x. « d, + 
(х,е) (je 3 το е) 
J” (x,€) = [jew : Є БЕ x mecs 


The direction finding problem is then modified to include the 
constraints in the sets Іс(х,е), I, (x,e), πε). апа J” (x,€). 


Therefore the direction finding problem becomes: 


Max: σ 
SUDI to: g - 7 Е(х) 15 < 0 
py (x)! « 0 περ x. E) 
x i`“ : Cr? 
al s < 0 Ее) 
—] — - 2 L =e 


ЕН 2 0, je) (x,e) 


2 - 
5; $0,7 (х,е) 
|5, | ШІ je (X) 


Nite solution to the direction finding problem is such 
that o < e, the problem is solved again with e/2. Thus we 
Stay away from the "corners" until we get near the optimal 
point. 
Іле бер Length Problem 
Once the direction = парсел found, the problem 
ID cu m hew far to move in that direction. If λε 


ЕАО Бе the largest step length A such that xÉ + ДЕЙ 


14 





ШЕП К, them the step length will be A, such that F(x* + А 55) 


F: 1£ A, а пеш еще 


problem is unbounded, and we can conclude that'the maximum of 


is the maximum of F(x* + ask) A <A 


F(x) is infinite. The maximum feasible step Ag is found by 
a simple search along the direction sk from the point х^. 
Then the optimal step length A, 1s determined by finding the 
k. К 
5) 5 


ο ο tien λε из used for the step length, 


А such that v F(x“ + zu сопсе Е(х) 1s concave, ТЕ 


v F(x“ 2 сащ se 
otherwise AL can be found by regula falsi [2]. 
4. Stopping Conditions 

ο ο αμ ο πι] terminate whenever the solution 
ЕО (пе direction finding problem yields c - 0, indicating 
that there are no usable feasible directions. If the solu- 
tion to the step length problem is λκ = ®, then the problem 
is unbounded and the procedure is terminated indicating an 
unbounded solution. In addition, if the step length or the 
Percent Improvement in the objective function becomes smaller 


than user provided limits, the problem is terminated, and 


the current solution is specified. 


C. THE INITIAL FEASIBLE POINT 
If an initial feasible point is not known, pick ап х! 


such that 5 ο. 


р < σης and solve the problem: 


Max: EX 
sub), to: T ге 
а; х - р: 21 з зу а. lel, 
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1 - 
where ΠΤ χι «Ὁ 
1 - i 
pce 
| it Е 
СИЕ оса. χὶ > Dp.. 
=i = і =1 = і 


An initial feasible point for this linear programming problem 
ES Хх”, 21 = 1. ТЕ the feasible region is non-empty, then the 
optimal solution of this problem will be z, = 0, x = x". Now 
PS si feasible with respect to the linear constraints, 
ШО ers also feasible with respect to the non-linear con- 
15, then use x" for the initial feasible point. How- 
“chet x for the initial feasible point. However if x" 15 
not feasible with respect to the non-linear constraints, 


solve the problem: 


Max: "2, 
Subj to: G; (x) -Р;2. < b: : ielo 
T 25 
а; х < b; ; iel 
dex Tc 
eae ee J 
22 > 0 
where: 
“ - tt 
р = 


G, (x'') - b: ΠΕ б. (х") > b:. 
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This non-linear programming problem can be solved by the 


method of feasible directions using X", z, = 1 as the initial 


Beugcuble point, If the feasible region of the original prob- 


lem is non-empty, then the optimal solution to this problem 


will be z, = 0, and x = x’. Now the original problem can be 


Q 


solved using x as the initial feasible point. 
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Pre Мо ИЧОМРОТЕК PROGRAM FEASBL 


This section will describe the computer oran aspects 
of the code, a complete set of user instructions is provided 
in section IV. Several programming languages were considered 
for the program including IBM 360/FORTRAN IV, SIMSCRIPT 11.5, 
and PL/1.  FORTRAN IV was chosen because its virtually uni- 
versal availability and widespread familiarity will allow 
easy installation at other computer centers, and facilitate 
user modifications in the future. 

FEASBL is structured as a straight line program, consist- 
ing of several sections; input, direction finding problem 
building and solving, step length determination, and output. 
In deciding the size of the problem that the program would 
be able to handle, the size of the direction finding linear 
programming problem was the o ne Баспай, То тпопте 
that round-off errors would not degrade the solution, the 
inverse of the basis matrix is maintained in double precision. 
FEASBL is capable of handling fifty variables, with fifty 
non-linear and fifty linear constraints thus requiring a 
basis matrix that is 152 by 152, which requires over 184K of 
core. Desiring to keep the total size below 300K if pos- 
sible, substantial effort was made to minimize the storage 
requirements throughout the program both by not storing any 
value that could be computed when needed, and using the same 


storage space for multiple purposes. The code now requires 
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about 290K of core space in addition to that required for 
user provided functions. 

The information required by the code includes the problem 
description, the initial starting point, and some program 
control variables. The problem description is provided by 
Pp ocirying themenyoctmve function, non-linear and Linear 
constraints, upper and lower bounds on the variables, and the 
requirements vectors. In addition, the first partial deriva- 
tives of the objective function and the non-linear constraint 
functions are required. A prototype code was developed which 
did not require derivatives, but which severely limited the 
form of the functions that could be used. It was felt, how- 
ever, that the greatly increased power and flexibility gained 
by not restricting the form of the functions more than offset 
F requirement to provide the derivatives. These functions 
are provided in the form of Fortran FUNCTION routines, thus 
FEASBL can handle any function the user can express as a 
Fortran routine. The remaining data is input Iin the Fortran 
NAMELIST format. This method was chosen because the unstruc- 
tured format minimizes keypunch errors. 

Each iteration begins by building the direction finding 
problem. The sets Тс (х,=), IT, (x,€), Т (х,е) ара Е ате 
determined, and then the direction finding linear programming 
problem is built. Rather then store the entire A matrix, 
only those elements which depend on the gradients of the ob- 
jective function and non-linear constraints in Ie (XE) referred 


to as the core are stored. The linear programming problem is 
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solved using the revised simplex algorithm [3], and the basis 
inverse is maintained through the product form. At each step, 
every non-basic variable is priced, and that variable with 
the most negative 2. š 5 ls selected for cme απ 
Sent of ties the lest-most variable is picked. The column 
vector а к for the entering variable is then generated, and 
used to compute y, - e ау. The departing vamdlable 15 j 
such that με. is minimum for all "E > 0. ІТ there 

are no negative 2; = S D then the optimal solution has 

been found for the direction finding problem. After an op- 
timal solution has been found, the basis matrix is created 
and inverted directly using the IBM subroutine DMINV, and 

the solution is checked for optimality. ТЕ the solution is 
not optimal, then the linear program is continued until an 
optimal solution is reached. When the optimal solution is 
achieved, the direction is available, and the program pro- 
ceeds to find the best step length. This is accomplished by 


k k 


first finding À the largest XA such that x + As? is fea- 


F? 
Able and then determining Лк < АЕ that maximizes F(xM+A, 5°). 
The quantity Ar is found by starting with A = 1 and increas- 


me А Бу а factor of 2 until aA is found such that 


NF 
xX * Àwp sk is not feasible. The interval between ХА) and 
the largest known feasible step Ar is reduced by checking 


À = (АЕ + An) ο. οἱ δι ιο until the interval is Tess 
then the minimum step length (STPMIN). However, if a feasi- 
pile step length A is found such that ER + 55) < F(x¥) , 


then it is assumed that the maximum point along the direction 
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lies between these two points. Once λε has been determined 


x is computed by the method Of regula falsi. Since F(x) is 


concave and is increasing at x" а ГОН НЕТ ІТ СЕ ко ο΄ the 
directional derivative (v F(x)" s) is positive at = and will 


be zero at the maximum along ss. If the directional deriva- 


mive 1s negative at x“ + Mas. then the maximum lies between 


the two points. A new À is then selected such that Л = 


(by - a\,) / (b - a) where: a = v F(x“ + 11569); b = v F(x“ 


1 


% ШЕТ; апа hy and A, bracket the optimal stop length. If 


2 


the directional derivative is zero at the new point the maxi- 
mum has been found, otherwise we have a new interval with 
positive and negative directional derivatives at the end- 
points, and another A is determined. This continues ипї11 
either the maximum is found, or the interval becomes less 


then the minimum specified. If the directional derivative 


is positive at xÉ + Aus. and NES Í Ans) > (хе), thempthe 


maximum along sk 1$ outside the feasible region, and Ху 15 


set equal to A However, if Е (хк + Aus) < F(x), then the 


F ° 
function is not concave, so a step hy between 0 and Ar 15 
found such that the directional derivative is negative, and 


MOS assumed that F(x) is concave between xÉ and хк B eres. 
and the maximum is sought as before. 
Once a usable feasible direction, and step length have 


been determined set ale = хк + MR and repeat the process 


= 
Starting at xÉ E. σος processus continued until either 
no usable feasible direction can be found or the problem 


becomes unbounded. If the step length or the percent 
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ENDrovement in the objective function become less then user 
supplied limits the program will terminate. The program will 


also stop after a user specified number of iterations. 


Да 





ТҮ. USING FEASBL 


This section will describe how to use FEASBL, including 
required input, possible errors, and the normal output. An 


example is given illustrating different features of the code. 


A. - INPUT 

The objective function, non-linear constraint functions, 
and their first partial derivatives are input as Fortran 
τον routines, The linear constraint matrix, require- 
ments vectors, upper and lower bounds on ze te ΠΙΕΙ ΕΞ. 
sible point, and any other input is read in as data at 
Execution. 

1. The FUNCTION Routines 

Four Fortran FUNCTION routines are required, these 

nEoutInes provide the objective function, its first partial 
derivatives, the non-linear constraints, and their first 
partial derivatives. Even if the problem has no non-linear 
constraints, all four FUNCTION routines must be provided. 


The four routines are: 


FUNCTION Definition 

FX( X ) The objective function. 

DELFX( NX, X ) The partial derivatives of 
the objective function. 

GIX( NGI, X ) The non-linear constraint 
functions, 


DELGIX(NX,NGI,X) The partial derivatives of 
the non-linear constraints. 
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The arguments of the FUNCTIONS statements are: 


Argument Definition 
X ШЫ vector x. 
NX Subscript of the partial 


derivative desired, i.e. 
IF (X) /9Xyy 


NGI The non-linear constraint 
number. 

Each routine must include a statement dimensioning X(e.g.; 
REAL*4 X(8), or DIMENSION X (15) ), at least one assignment 
statement with the name of the FUNCTION on the left of the 
а (ер FX = ХТ), or DELGIX = 4.0 * X(2)**3 ), 
and at least one RETURN statement. No variation in the names 
of the FUNCTIONs, nor in the order of the arguments is ai- 
Towed. The following example illustrates the use of the 
FUNCTION routines. 
2 2 
1 


Мах: IES = OR OX, тока я X» 


Subj to: 4х2 + х < 16 
non-linear constraints 
(4 - X1) + τ < 16 


SX Хх < 15 | 


mear constraints 
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For this problem the following FUNCTION routines are required: 


a. ¡Me "objective function 


FUNCTION FX( X ) 

REAL*4 X(2) 

[ = б.0 * X(1) + 8.0 * X(2) = X(D 2 - xO: 
RETURN 

END 


b. The partial derivatives of the objective function 


101 


102 


FUNCTION DELFX( NX, X ) 
REAL*4 X(2) 

GO TO ( 101, 102 ), NX 
DELFX = 6.0 - 2.0 * X(1) 
RETURN 

DELFX = 8.0 - 2.0 * X(2) 
RETURN 

END 


c. The non-linear constraint functions 


100 


200 


FUNCTION GIX( NGI, X ) 

REAL*4 X(2) 

ало 100, 200 ет 

ВХ = 4.0 * Ха + Хо] 82 
RETURN 

X ы AA) кк 
RETURN 

END 


d. The partial derivatives of the non-linear con- 
straints 


FUNCTION DELGIX( NX, NGI, X ) 
REAL*4 X(2) 
ЕШ слова 2000 NGI 


ІШІНЕО TO ( 101. 1025 М 
NOUBDERGIX = 8.0 * X(1) 


RETURN 


ША ШЕТСІЗ = 2.0 * X(2) 
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RETURN 

200 GO TO ( 201, 202 ), NX 

ПОЛЕ СЕО πι] = 8.0 
ВЕТОВМ 

202 DELGIX = 2.0 * X(2) 
RETURN 
END 


The significant feature to note is that when a FUNCTION is 
called by FEASBL specifying a point X, and when necessary 
the partial derivative desired NX, and the constraint number 
NGI, the appropriate value is returned. 

2. Data Input 

The remainder of the data needed is read in from 

data cards at execution in NAMELIST format. This format was 
chosen to minimize Keypunch errors, anda provide flexibility 


for the user. The following variables are required: 


Name Description 

IXMAX Number of variables. 

IGCNT Number of non-linear constraints. 
ILCNT Number of linear constraints. 


The following variables are optional, and if not provided by 


ene user will take on the default values listed. 


Name Default Description 

GL(I,J) 0.0 Linear constraint matrix element 

XMIN(I) 0:0 Lower bound on E 

XMAX (I) :10/? Upper bound оп х, 

BC(I) 0.0 The bi for the поп- 1 саг con- 
stralımts. 
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Name Default Description 


BL(I) 0.0 The b: for the linear constraints, 

X(I) 0.0 The initial feasible point. 

EPSLON 9.1 The initial anti-jamming constant. 

STPMIN 10% The minimum step length. 

PERCENT D The minimum objective function 
1mprovement. 

ITMAX 25 The maximum number of iterations. 


The first data card begins with an ampersand (8) in column two 
followed by the word INPUT with no space ( &INPUT ). Тһе 
values are then listed as "variable name = FORTRAN constant" 
1.е. ІХМАХ - 2, GL(1,1) = 3.0, the last data item is followed 
by an ampersand and the word "END" ( &END ). If more then 

ИС сата is required, each сага must begin in column 2 or 
more, and the ¿END will appear on the last card only. For 


the above example the data cards would look like: 


GINPUT IXMAX = 2, IGCNT=2, ILCNT=2, GL(1,1)=3.0, GL(1,2)=5.0, 
БІН  -17)-2:0, 62(2,2)-12:0; ХМАХ(1)-1.85, ХМІМ(2)-2.0, 
В!) -16.0, ВС(2)-16.0, В.(1)-15.0, В.(2)-30.0, ITMAX=10, 
ВЕ 0. 01, ῬΙΡΜΙΝΞΙ.0Ε-04, X(1)=1.0, X(2)=2.0, GEND 
Here.the order is unimportant, and all the variables need not 
be listed. In this example the default values for XMIN(1), 
XMAX (2), and PERCNT will be used since they are not listed 
on the data cards. 
5. Data Limits 
Input errors will result if the variable inputs arë 


outside the following limits: 
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1 < ТХМАХ < 50 


c 
^ 


IGCNT « 50 


© 
A 


ILCNT < 50 
XMIN(I) « XMAX(I) 


To prevent stopping prematurely, EPSLON should be significantly 
larger than STPMIN. 

After all the data is input, the initial point is 
puecked for feasibility. If the point is not feasible the 
program is terminated. 

4. The Program CHECK 

FEASBL requires 290K of core space in addition to 
that required for the FUNCTION routines. Most problems will 
require less then 320K and less then 4 minutes to run. As 
an aid for checking the input data and FUNCTION routines, the 
program CHECK can be used. CHECK requires exactly the same 
input as FEASBL, and performs the same input error checks. 
In addition CHECK will determine the feasibility of the ini- 
τὴ Ρο απ and print out the value of the objective function, 
all the constraints, and the value of all the partial deriva- 
ето the Objective function and non-linear constraints. 
By using CHECK the user can quickly verify the input, since 
CHECK can be run normally as a class B or C job, before 
Using FEASBL, which is a class J job. 

5. system Cards 
The JCL cards required to use FEASBL or CHECK depend 


on the way that the FUNCTION routines are input. The routines 
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can be input as either Fortran source decks and compiled as 
ШЕР ОЕ the Job, or the routines can be precompiled and in- 


put as object decks. An example of each method follows. 
а. FEASBL with Source Decks 


// standard OS JOB card ‚ТТМЕ= (05,59) 
// EXEC FORTCLG,REGION.GO=300K 
//FORT.SYSIN DD * 
(FUNCTION source decks) 

/* 
//LINK.NLPLIB DD DSNAME=F0600.NLPLIB,UNIT=2321, 
//  VOL=SER=CEL001,DISP=SHR,LABEL=(,,,IN) 
ЕРМЕ SYSIN DD * 

INCLUDE NLPLIB(FEASBL,DARRAY) 

/* 
"SO SYSIN DD * 

GINPUT IXMAX-2, IGCNT=2, ... , GEND 

/* | 


29 





b. FEASBL with Object Decks 


7 standard 0S JOB card ТЕ ου 
// EXEC LGO,REGION.GO=300K 

//LINK.NLPLIB DD DSNAME=F0600.NLPLIB,UNIT=2321, 
// | VOL-SER-CEL001,DISP-SHR,LABEL-(, , , IN) 
ГАМК. УБТ DD * 

INCLUDE NLPLIB(FEASBL,DARRAY) 
(FUNCTION object decks) 

/* 
ГОО БИЛМ DD * 

&INPUT IXMAX=2, IGCNT=2, ... , X(2)=2.0, &END 
/* 
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с. CHECK with Source Decks 


// standard OS JOB card 
// EXEC FORTCLG,REGION.GO=100K 
//FORT.SYSIN DD * 

(FUNCTION source decks) 
/* 
//LINK.NLPLIB DD DSNAME=F0600.NLPLIB,UNIT=2321, 
// VOL=SER=CEL001,DISP=SHR,LABEL=(,,,IN) 
ЖЕТІК БҮСІМ рр 5 

ВОВЕ КЕРЕТВ (| СНЕСК) 
/% 
MAGO. SYSIN DD * 
&INPUT IXMAX=2, IGCNT=2, ... ,X(2)=2.0, SEND 

/* 
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da CHECK wren Ob Cece Decks 


ЙА standard OS JOB card 
// EXEC LGO,REGION.GO=100K 
//LINK.NLPLIB DD DSNAME=F0600.NLPLIB,UNIT=2321, 
// NOL-SER-CEL001,DISP-SHR,LABEL-(, , , IN) 
//LINK.SYSIN DD * 
ТСИ ЕРЕ Е Е) 
(FUNCTION objective decks) 

/* 
//GO.SYSIN DD * 

&INPUT IXMAX=2, 


/* 


IGCNT=2, , X(2)=2.0, &END 


A set of object decks can be produced either atter 


checking the FUNCTION routines by executing a FORTCD job, or 


in conjunction with the program CHECK by replacing // EXEC 


FORTCLG, 


REGION.GO=100K with // EXEC FORTCLGP,PARM.FORT- 


'"DECK',REGION.GO-100K. 


Br 


The 


SAMPLE PROGRAM 


FUNCTION routines, input cards, and system cards 


required to use FEASBL with source decks are shown for the 


Problem: 


| В Жа 2 
Мах: ρω.» OX, + 8x, X1 x) 
Subj to: 4х2 + х2 < 16 

1 2 - 
2 2 
(4 - X1) + x, < 16 
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5X. + 5X £5 


tA 


ZA Лх 50 


= 
wN 
ΙΛ 


A sample of the output of both CHECK and FEASBL for this 

Euoblem rs provided. The next two pages are the cards re- 
quired for executing FEASBL with the sample problem. This 
is followed by the output of CHECK and then FEASBL for the 


sample problem. 
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REAL*4 X(2) 
THIS 15 THETOBJECTIVE FUNCTION 
= en BRUTTO О) IE 2 =< X(2)]*%*2 


ας ας ας ἃς о a fe О Хо О ОО ЖЖ ake ok 


Х 
E THE FIRST PARTIAL DERIVATIVES 


) 
1, 102 ), NX 

> THE PARTIAL WITH RESPECT TO X(1). 
Ξ 6.0 - 2.0 * X(1) 


X 
R 
THIS IS THE PARTIAL WITH RESPECT TI Х(2) 
2 = 8.0 - 2.0 ж Х(2) 


FUNCTION GIX( NGI, X ) 
THIS FUNCTION PROVIDES THE NON-LINEAR CONSTRAINTS. 
REAL*4 X(2) | | 
GO TO THE PROPER CONSTRAINT. 
GO TO ( 100, 200 ), NGI 
THIS IS NON-LINEAR CONSTRAINT NUMBER 1. 
X il * X(1)**2 ж Х(2)хж2 


THIS IS NON-LINEAR CONSTRAINT NUMBER 2. 
( 4.0 - X(1) )жж2 + х(2)»жж2 


| 


GIX = 
RETURN 
END 
жек к со sk oie ood tool эсс Же с ЖЖ ккк Ж Жк Жок кок жокк Жжжж ж 
FUNCTION DELGIX( NX, NGI; X 


πο ο νετ PROVIDES ΤΗΕ 
OF THE NCN-LINEAR CONSTFAT 


) 
FIRST PARTIAL DERIVATIVES 
NT FUNCTIONS. 
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100 
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200 


ООО ООО 


201 


ООО 


202 


REAL*4 Х(2) | 
GO TO THE CONSTRAINT SPECIFIED. 
GO TO { 100. 200 ),NGI 
THIS IS CONSTRAINT NUMBER ONE, 
GO TO ( 101, 102 ), NX 
THIS IS THE PARTIAL WITH RESPECT 
πῇ = 8.0 * X(1) 
THIS IS THE PARTIAL WITH RESPECT 
DELGIX = 2.0 х Х(2) 
RETURN 


WHICH 


DELG 
RETU 


THIS IS CONSTRAINT NUMBER TWO, 
СО ТС { 201, 202 ), NX 
THIS IS THE PARTIAL WITH RESPECT 


IX = 2.0 * Х(1) - 8.0 
RN 


IS THE PARTIAL 
2.0 * X(2) 
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TO X(1). 


TO X(2) 


PARTIAL IS WANTED 


TO X(1) 
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SAMPLE OUTPUT FROM CHECK 


NON-LINEAR PROGRAMMING PROBLEM 


INITIAL SET UP: 
NUMBER OF VARIABLES: 
NUMBER OF NON-LINEAR CONSTRAINTS: 2 
NUMBER OF LINEAR CONSTRAINTS: 
THE LINEAR CONSTRAINT MATRIX IS: 
C1, 1}: 3.00000 ( 1, 2); 5.00000 


( 2, 1); 2.00000 ( 2, 2); 12.0000 


THE BCUNOS ON THE VARIABLES ARE: 


0.0 < ХЕ 1) < 1.85000 ; 2.00000 «πι г) < О. 723701Е Το 
NON-LINEAR CONSTRAINTS REQUIREMENTS VECTOR, BC(I): 
( 1); 16.0000 ( 2); 16.0000 ( 
LINEAR CCNSTRAINTS REQUIREMENTS VECTOR, ВЕЕТ): 
(1); 15.0000 (2); 30.0000 ( 
THE INITIAL FEASIBLE STARTING POINT IS: 

Xt 1); 1.00000 ‚Xt 2); 2.00000 : Xt 

ANTI JAMMING CCNSTANT: 0.160000c-01i 
MINIMUM STEP LENGTH: 0. 100000E- 03 
MINIMUM IMPROVEMENT: 0.100000E-01 
MAXIMUM 4 ITERAYICNS: 10 


INPUT CCMPLETE. INITIAL VALUE OF OBJECTIVE FUNCTION: 1. 700000Е 01 


VALUE CF NON-LINEAR CONSTRAINT NUMBER, 1, IS: В.0000000Е 90 

VALUE Cf NON-LINEAR CONSTRAINT NUMBER, 2, 15: 1. 3000000Е 01 

THE VALUE OF LINEAR CCNSTRAINT NUMBER І, 15: 1.300900E Oi 

ΤΗΕ VALUE OF LINEAR CCNSTRAINT NUMBER 2, 15: 2.6090000E Ol 

THE BCUNOS, ANO ALL THE NCN-LINEAR ANO LINEAR CCNSTRAINTS HAVE ВЕЕМ СНЕСКЕО. 
NEXT THE VALUES OF THE PARTIAL OERIVATIVES OF THE DBJcCTIVE FUNCTION ANO THE CONSTRAINTS ARE LISTEO. 
THE VALUE OF OELFXE 1, X) IS: 4.000000E 00 

THE VALUE CF OELFX( 2, X) IS: 4.000000Е 00 

THE VALUE OF OELGIX( 1, 1, XI IS: B.000000E 00 

THE VALUE OF OELGIX( 2, 1, X) IS: 4.090000Е 00 

THE УА(МЕ ОҒ ОЕЦСІХІ 1, 2, Х) 15: -6.000000Е 00 

THE VALUE OF OELGIX( 2, 2, XI IS: «.000000Е 00 
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SAMPLE OUTPUT FROM FEASBL 


NON-LINEAR RROGRAMMING PROBLEM 


INITIAL SET UP: 
NUMBER DF VARIABLES: 
NUMBER OF NON-LINEAR CONSTRAINTS: 2 
NUMBER OF LINEAR CCNSTRAINTS: 
THE LINEAR CCNSTRAINT MATRIX IS: 
μιν: 3.00000 i25 5.00000 


( 2, 1); 2.00000 ( 2, 2); 12.0090 


THE BSUNOS CN THE VARIABLES ARE: 


0.0 < X( 1) < 1.85000 ; 2.00000 < X( 2) < 
NON-LINEAR CONSTRAINTS REQUIREMENTS VECTOR, BCCI): 
(1); 16.0000 ка); 16.0000 ( 
LINEAR CONSTRAINTS REQUIKEMENTS VECTOR, BL (I): 
2 117; 15.0000 ( 2); 30.0000 ( 
THE INITIAL FEASIBLE STARTING ROINT IS: 

Xt 11; 1.00000 »X( 2); 2.00000 ; Xt 

ANTI JAMMTNG CONSTANT: 0.199009=-01 
MINIMUM STER LENGTH: 0.100000£-03 
MINIMUM IMPROVEMENT: 0.100000£-01 
MAXIMUM # ITERATICNS: 10 


INPUT COMPLETE. INITIAL VALUE OF OBJECTIVE FUNCTION: 1.700000€ Ol 


PRCBLEM SOLUTI^N 


ТЕКАТ NUMBER: І 
pe Fix) = E “1 867500€ 01; IMPROVEMENT THIS ITERATION: 11.0294 RERCENT. 


THE ROINT X( ) IS: 


Xt I) 1.250000 2 Х( 2), 2.250000 $ Х( 
NR 2112 Ol: IMPROVEMENT THIS ITERATION: 1.8393 PFRCENT. 
THE POINT X( ) IS: 
X( 1), 1.666625 > XO 2), 2.000024 $ Х( 
CRTIMAL SOLUTION: 
THE ORTIMAL VALUE CF THE CBJECTIVE FUNCTION IS: 1.922217€ Οἱ 
TH£ POINT X( ) IS; 
X( 1), 1.666625 ; X( 2), 2.00002% $ Х( 


oy 
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V. EVALUATION 


To validate the code, and evaluate its performance, 
several non-linear problems were solved with FEASBL. These 
problems were selected to represent the range of problems 
FEASBL was designed to solve. The problems included uncon- 
strained optimization of a concave objective function, con- 
strained optimization of concave objective functions some 
with the unconstrained maxima interior and some exterior 
to the feasible regions. An unbounded problem was also 
solved. 

In each case the optimum solution was obtained in less 
пишеш five сесопйе ос ері үүле, Іневе ргор1ете Mere л тшс ΠΠ 
in order to facilitate manual verification of the solutions: 
Larger problems will take correspondingly more time, but the 
code appears to be quite efficient. In most problems the 
optimal solution was obtained after two iterations. Іп опе 
case, the solution took seven iterations. 

улг оп εμας πε value orf STPMIN is critical to the 
performance of FEASBL. If STPMIN is set too small, when the 
interval between Ar and АЕ ΤΙ. ἘΠῈ round ott 
OL may exceed the value of STPMIN. This results in the 
code looping indefinitely trying to reduce the interval beyond 
machine accuracy leading to a program completion due to ex- 
cess time. If this should occur the problem should be re- 
Started at the last point with a amines value του ο Е 


With step lengths on the order of 1.0 the interval cannot be 
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determined to an accuracy greater then 10 9. With larger 


step lengths this value will increase. An initial value 


2 Or 107? is recommended for STPMIN. This value 


of 10. 
cam be reduced in subsequent runs of FEASBL if a more pre- 
cise solution is desired. 

The initial value of EPSUON is not as critical as STPMIN. 
EPSLON should be greater than STPMIN but generally not greater 
1 


than 10° Experience with FEASBL indicates that values on 


1 or 1072 


the order of 10 perform well. If a large value is 
used, on the order of 1.0 or greater, and an optimal point 
is found, this should be verified using a smaller value for 
EPSLON, as the code may stop early with large values of 
EBSLON. 

Performance was uneffected by variations in PERCNT. Large 
values may cause the program to terminate sooner. Values less 
then 1079 are equivalent to zero because of inherent machine 
accuracy. 

An attempt was made to solve a problem in which all the 
constraints were equality constraints, thus the feasible 
region had no strictly interior points. The program would 
terminate immediately without any significant movement. This 
ПОПЕ was anticipated since the underlying theory assumes 


that the feasible region contains at least one strictly in- 


ιτ pont to the set of non-linear Constraints [1]. 
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